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Exclusive weak radiative B meson decays are studied using the covariant oscillator quark 
model. The branching ratios for the processes B" K*°-i and B± ^ K'^-y have been 
estimated. These are in reasonable agreement with the available experimental data. The 
^ calculation has been extended to the CKM-suppressed decay processes B~ ~* p^j and 

in . ^ if *''7. 
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I §1. Introduction 

o 

. The study of the weak radiative decay B K*^ as a test of Standard Model 

(SM) has attracted considerable attention since the CLEO experimental gave the 
Q_i' preliminary determination of the exclusive branching ratio Br(B — > K*^) = (4.0 it 

1.7 ±0.8) X 10~^. The weak radiative decays of B mesons (which proceed through a 
QJ \ flavor changing neutral current, absent at the tree level in the SM) are remarkable for 

' several reasons. The B — > K*'y decay arises from the quark level process 6 — > 57 via 

penguin-type diagrams at the one loop level. Hence it is not only a significant test of 
^ I Standard Model fiavor-changing neutral current dynamics but also is sensitive to new 

^ ' physics appearing through virtual particles such as the top quark and W boson in 

the internal loop. The study of this process provides valuable information concerning 
the Cabibbo-Kobayashi-Maskawa (CKM) parameters, Vtd, Vts and Vtb- Furthermore, 
additional contributions in loop stemming from new bosons and fermions present in 
most of the extensions of the SM, suggests the possibility that this process provides 
a window to new physics. 

The theoretical analysis of the decay process B — > K*^ requires long-distance 
QCD contributions which cannot be determined perturbatively. It is also not straight- 
forward to calculate the exclusive decays by first principles of QCD, due to compli- 
cations inherent in nonperturbative QCD. Therefore, one must resort to some phe- 
nomenological models to obtain reliable results. The heavy quark effective theory 
(HQET)lP is expected to be useful in this regard in so far as the b quark is con- 
cerned. However the s quark in the final hadron can neither be considered heavy 
enough to enable the use of HQET nor sufficiently light to permit the exploitation 
of the chiral perturbation theory in an unambiguous manner. There are several 
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methods available in the literature to stucLylhe exclusive B —^K*j decay process. 
Some of them include theJX^D siLni_t;ule,a^'Q^'l3) lattice QCD,cP nonrelativistic and 
relativistic quark models. The HQExlli has also been appHed to this 

decay process even though the s-quark mass is certainly not heavy enough, contrary 
to the requirement of HQET. 

In this paper we study the rare radiative decay process B — > K*'y using the 
covariant oscillator quark model (COQM).!!^ One of the most important motives 
for this model is to describe covariantly the center of mass motion of hadrons, while 
preserving the considerable success of the non-relativistic quark model regarding 
the static properties of hadrons. A keystone in COQM for doing this is treating 
directly the squared masses of hadrons in contrast to the mass itself, as done in 
conventional approaches. This makes Jiie covariant treatment simple. The COQM 
has been applied to various problems EJ^ with satisfactory results. Recently Ishida 
et al.Ej''ll3) have studied the weak decays of heavy hadrons using this model and 
derived the same relations of weak form factors for heavy-to-heavy transition as done 
in HQET.B^ In addition, our model is also applicable to heavy-to-light transitions. 
As a consequence, this model does incorporate the features of heavy quark symmetry 
and can be used to compute the form factors for heavy-to-light transitions as well, 
which is beyond the scope of HQET. Actually, in previous papers we made analyses 
of the spectra of exclusive semi-leptonic EJ) decays of i?-mesons and analyses of non- 
leptonic decays of B mesons and of hadronic weak decays of Ab baryonsEfP using 
this line of reasoning, leading to encouraging results. Keeping this success in mind, 
we extend its application to weak radiative decays of B mesons. 

The paper is organized as follows. In §2 we present the expressions for the decay 
widths for the weak radiative decays of B mesons. In §3 we present a brief description 
of the covariant oscillator quark model. Using this model we have evaluated the 
required form factors. Section 4 contains our results and discussion. 



The general amplitude of weak radiative decay with one real photon emission is 
given by 



where r] and e are the polarization vectors of the photon and the vector meson P* , 
respectively. The first (second) term on the right-hand side of Eq. (2-1) is a parity 
conserving (violating) term, and because of the real photon we have g^=0. The 
decay width implied by this amplitude is given as 



§2. Methodology 



+ t U^{Ml -Ml.)-{p + k)^{e ■ q)] ^(g^) 



(2-1) 
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1 (M| -M^p^f 



|/i(0)p + 4|/2(0)p 



(2-2) 
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Now, the exclusive decay B K*^ is expected to be described well by the quark 
level process 6 — > 57 if the confinement effect is properly taken into account. Accord- 
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ingly, in the Standard Model, B decays are described by the effective Hamiltonian 
obtained by integrating out the top quark and W boson fields given as 113) 



AGf 



Vtb Vtl E ^^(/^) , 



(2-3) 



where the Ci{fj,) are the Wilson coefficients, arise from the renormalization group 
equation to provide the scaling down of the subtraction point appropriate to the 
problem i.e., /x ~ ruf,. The set {Oj} is a complete set of renormalized, dimension-six 
operators which govern the b ^ s transition. They consist of two current operators 
Oi and O2 and four strong penguin operators O^-Oq, which determine the nonlep- 
tonic decays, and the electromagnetic dipole operator O7 and the chromomagnetic 
dipole operator Og, which are responsible for rare B decays, i.e., 6^5 + 7 and 
b ^ s + g, respectively. Out of the eight operators Oi, the one that contributes to 
B K*-f is 



O7 



327r- 



:Ff,^ [mbs a'"' (1 + 75) & + a^"" (1 - 75) b] , 



(2-4) 



where Ffj^i, is the electromagnetic field strength tensor. The relevant Wilson coeffi- 
cient is given bytH^ 



C7(/x) 



as{mh) J I 

/ N X 2/23 



CriMw) - -Cs{Mw) 



with 



and 



CriMw) 



as{Mw). 
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Qs("T-b) 

a, (Ml 



w 



19/23 



(X - 1)3 



CsiMw) 



\J2 
2 



+ 



[X - 1)4 

3xlnx 



{x - 1)3 [x - 1) 



(2-5) 



(2-6) 



(2-7) 



where x = m^/M^. From the fact that ^ ms, only the term involving mf, in 
the operator O7 need be retained. Thus the matrix element of interest becomes 



(i^*7|07|S) 



le 
16^ 



q^^T], ^KXso^"" (1 + 75) h\B) 



where is the four momentum of the photon and t/^ is its polarization vector. It 
should be noted that here the matrix element of a tensor current between hadronic 
states for which not much information is available is involved. However, in the frame- 
work of HQET the associated heavy quark spin symmetry enables one to express 
the matrix element of the tensor operator in terms of the vector and axial vector 
form factors that also occur in semileptonic decays and may be estimated in different 
phenomenological models. 
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In the static limit pi a heavy b quark we may use the equation of motion 796 = b 
to derive the relationsE2P 

{K*\s i aoi (1 + 75) b\B) = {K*\s 7, (1 - 75) b\B) . (2-9) 

As a result, the form factors /i and /2 in Eq. (2-1) can be related to the vector and 
axial vector form factors V and Ai appearing in the matrix element of the RHS of 
Eq. defined by§) 

{K*{k)\sj,b\B{p)) = ,^qyP v{q') , (2-10) 



{K*ik)\sj^j5b\Bip)) = (Mb + MK*)e^ Aiiq^) 



e ■ q 



Mb + Mk' 



2'-^q^MK* [A,{q') - Ao{q')] , 



with q = p — k. Here A^{q^) is simply an abbreviation for 

2^ Mb + Mk' . , 2. Mb - Mk* . , 2, 

^3{Q ) = ^1(9 ) TTTJ ^2(g ) , 



(2-111 



(2-12) 



and in order to cancel the singularity at = 0, we must have ^3(g^ = 0) = Ao{q'^ = 
0). Thus with Eqs. (2-1) and (^-( ^) , we obtain the following relations for 
B -> K* transition at q^ = 0:© 



/i(0) 



Gf e 
722^ 



C^iiS) VtbV:, F(0) 



and 



where 



/2(0) = -/i(0), 



Mb- Mk' . ^ Mb + Mk' . , . 
F(0) = ^ ^ y(0) + ^T. ^i(O)- 



(2-13) 



(2-14) 



(2-15) 



1Mb "^Mb 

Substituting the above values of /i and /2 into Eq. (|2-2D , we obtain the decay width 
for B K*'j as 



1287r4A/5 



(Mb + Mk')Ai{0) + (Mb - M^O^W 



(2-16) 



Similarly, the decay width for the CKM suppressed FCNC radiative transition 
b — > d'j, which is responsible for B~ p~'j and Bs K*^, is obtained from 
Eq. (2-16) with Mb and AIx* replaced by the corresponding initial and final meson 
masses. The relevant CKM factor in this case is |VtbV((ip. 
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Now to evaluate the form factors ^i(O) and V{Q), we use the covariant oscillator 
quark model, which is presented in the next section. 

Here it should be noted that the relations ( 2-13 ) and ( 2-14 ) are also derivable Ej) 
directly from Eq. ( |2-8| ) in COQM, by using the covariant spin wave functions, that 
the Bargmann-Wigner spinor functions. Accordingly, the expression of -F(O) in 



IS 



terms of a space-time wave function (Eq. (3-10) given later) is also derived directly. 



§3. Model framework and the hadronic form factors 



The general treatment of COQM may be called the "boosted -LS'-coupling scheme," 
and the wave-functions, being tensors in C/(4) x 0(3, l)-space, reduce to those in 
5'C/(2)spin X 0(3)orbit-space in the nonrelativistic quark model in the hadron rest 
frame. The spinor and space-time portion of the wave functions satisfy separately 
the respective covariant equations, the Bargmann-Wigner (BW) equation for the 
former, and the covariant oscillator equation for the latter. The form of the meson 
wave function has been determined completely through the analysis of mass spectra. 

In COQM the meson states are described by bi-local fields ^^(xi^, X2/x)) where 
xi^{x2^j) is the space time coordinate of the constituent quark (antiquark), A = 
(a, a) {B = (6, (5)), describing its flavor and covariant spinor. Here we write only the 
positive frequency part of the relevant ground state fields: 

^2m) = e'''-'' U{P)i fabix^; P) , (3-1) 

where U and / are the covariant spinor and internal space-time wave functions, 
respectively, satisfying the Bargmann-Wigner and oscillator wave equations. The 
quantity (X^) is the relative (CM) coordinate, = xi^ — {Xfj, = mixi^ + 
m2X2fj,) / {mi + 1712), with rrii the quark masses). The function U is given by 

UiP) = ^ [{-l5Psiv) + i7^.V^iv))il + zv . 7)] , (3-2) 

where Ps(V^) represents the pseudoscalar (vector) meson field, and = P^/M {P^{M) 
is the four momentum (mass) of the meson). The function U, being represented by 
the direct product of a quark and antiquark Dirac spinor with the meson velocity, 
is reduced to the non-relativistic Pauli-spin function in the meson rest frame. The 
function / is given by 

/(-.;^) = ^xp(-f + . (3.3) 

The value of the parameter j3 is determined from the mass spectrail' as (3p = 0.14, 
Pk* = 0.142, pB = 0.151 and 0.160 (in units of GeV^ ). 

The effective action for weak interactions of mesons with 1^-bosons is given by 



Sw = I d'^Xid'^X2{^F,P'{xi,X2)i'^^,{l+-1^)^I,p{xi,X2))W^^q{xi] 



(3-4) 
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where we have denoted the interacting (spectator) quarks as 1 (2) (the KM matrix 
elements and the couphng constants are omitted). This is obtained from consid- 
eration of Lorentz covariance, assuming a quark current with the standard V — A 
form. In Eq. (p-4|), ^i^p {^f,p') denotes the initial (final) meson with definite four 
momentum {P'p)-, and is the momentum of VF-boson. The function ^ is the 
Pauli-conjugate of <P defined by ^ = — 74(^^74, and ( ) represents the trace of 
Dirac spinor indices. Our relevant effective currents J^{X)pi^p are obtainedll^ by 
identifying the above action with 

Sw = j d^XJ^{X)p,,p W^{X)g, (3-5) 



+Vx{v')Ps{v){e^Xaf3 < Vf3 - 6x^^{w + l)-vx v'^] , (3-6) 

where M(M') denotes the physical mass of the initial (final) meson and w = —v ■ 
v'. The quantity I'^^{w), the overlapping of the initial and final wave functions, 
represents the universal form factor /unctionH It describes the confinement effects 
of quarks, and is given by 



4(3(3' 



: exp(-G(u;)) ; C{w) = (/? - P'f + 4(3 (3' w"^ 



(3 + (3' ^/C(^\ 



(3-7) 



and 



G{w) 



2C{w) 
+2 1/3' 



M 



M 



i\ 2 



M 
Ml 



ML 



[w 



MM' 
2 w 

MsM'^ 



1) 



(3-8) 



where Ms{M'g) represents the sum of the constituent quark masses of the initial 
(final) meson, and m„ is the spectator quark mass. 



Comparing our effective current (3-6) with Eqs. ( p-10 ) and (2-11), we obtain 
the relations between the invariant foriiLfactors V and Ai with the form factor 
function liw) for B K* transitions asli^ 



y(o) = yii(o) = F(o), 



m 



Mb + Mk* 



(3-9) 
(3-10) 



and similarly for B 



and Bs K* transitions with I'^^{'w) replacing P^{w). 



*' In this paper we apply the pure-confining approximation, neglecting the effect of the one- 
gluon-exchange potential. This approximation is expected to be good for the heavy/light-quark 
meson systems. 

**■* Here, note the remark given at the end of §2. 
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§4. Results and conclusion 

To estimate the branching ratios for weak radiative decays of B mesons, we 
use the following values for various quantities. The quark masses (in GeV) are 
taken as ruu = rud = 0.4, rus = 0.55 and rrib = 5. The particle masses and their 
lifetimes are taken from Ref . ^) . The relevant CKM parameters required for all the 
processes considered here are taken from Ref. ^) as Vtd = 0.0085, Vts = 0.0385 and 
Vtb = 0.99925. The renormalized Wilson coefficient Caijnb) used in the estimation 
of branching ratios is taken to be |C7(m;,)|=0. 311477. B^) With the above values we 
first evaluate the form factor F{0) using Eqs. ( |2-15| ) and (|3-7|)-(3-9). Our results 
for the form factor values given in Table I compare well with the recent calculations 
within the framework of the relativistic quark model, 1121' light cone QCD sum rule,cl^ 
and hybrid sum rule. EP With the calculated values of the form factors, the branching 
ratios for different channels are estimated using Eq. (2T6) as well as the mean 
life values of the appropriate decaying meson, which are given in Table II. The 
branching ratios Br{B^ — > K*^'j and Br{B^ —i- K*^j) are found to be in very 
good agreement with the available data. Here it may be worthwhile to note that 
the relevant process is relativistic, and the form factor function plays a significant 
role, such that / = 0.235. Since there are no data as yet available for the branching 
ratios in the case B ^ p-y and Bg — > K*'y, the model predictions must be compared 
with other theoretical predictions available in the literature. In fact there are few 
theoretical attempts made so far in this sector. Recently, SingercB* predicted that 
in the relativistic quark model, 

Br{Bs ^ K*-f) = {IA±0.6) X IQ-^, (4-1) 

and in the relativistic independent quark model, Barik et al.0) obtained 

Br{B pj) = 1.24 X 10"^, 

Br{Bs K*-f) = 9.28 x 10"^ . (4-2) 
Finally, the hadronization ratio R defined as r{B K*^)/r{h — > 37), where 

nb - «7) = ^ \VtsVtb? ml , (4-3) 

is found to be 0.12, which agrees well with the experimentally observed value R = 
0.19 ± 0.09.13) 

In view of the consistency of our predictions, obtained with no free parameters, 
with the large number of theoretical predictions as well as the experimental data, 
it may be concluded that the framework of COQM provides a suitable scheme to 
estimate the confined effect of quarks, and so far as the relevant flavor changing 
neutral current is concerned, the Standard Model is valid. 
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Table I. Prediction for the rare radiative decay form factor F{0) in comparison with various theo- 
retical predictions. 



Decay process 


Present work 


Ref. 


1C|) 


Ref. |4 




Ref. |5|) 




0.328 














0.329 


0.32 ±0.03 


0.32 ± 0.05 


0.308 ± 0.013 ± 0.036 ± 0.006 


B~ p"7 


0.295 


0.26 ± 0.03 


0.24 ± 0.04 


0.27 ±0.011 ±0.032 




0.250 


0.23 ± 0.02 


0.20 ±0.04 





Table II. Prediction for the branching ratios of the exclusive rare decays of B mesons along with 
their experimental values. 



Decay process 


Br ratio 
This work 


Br ratio 
Expt. p4|) 




3.96 xlO"^ 


(4.0 ±1.9) X lO"'^ 




4.168 xlO"^ 


(5.7 ±3.3) X lO"'^ 


B~ p"7 


1.68 X 10"'' 




Bs ^ /•i'*»7 


1.158 X lO"'* 
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